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One of the most debated issues related to high-T c super- 
conductivity is the symmetry of the Cooper pair or the gap 
function. In this report, we present numerical results regard- 
ing the gap function in strongly correlated electron systems 
using t — J and Hubbard models in one and two dimensions. 
To this end, we use exact diagonalization to study the ground 
states of 8- and 16-site clusters consisting of single or coupled 
layers. We calculate a reduced two-particle density matrix in 
momentum space which is a measure of the gap function. We 
then analyze the eigenvectors of this density matrix, which 
display the possible Cooper pair symmetries. The eigenvec- 
tor corresponding to the largest eigenvalue indicates a vanish- 
ing gap on the Fermi surface (which is in favour of odd-gap 
pairing) although d x 2_ y 2 symmetry is seen to be a very close 
contestant in many of the cases. 



I. INTRODUCTION 

In the theory of superconductivity, the quantity which 
characterizes the nature of pairing is the gap function. 
More precisely, it is the amplitude for pairing 
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for a Cooper pair of electrons in states (k T/,— k J.) in 
a singlet spin configuration. This quantity is closely re- 
lated to the energy gap in the single-particle excitation 
spectrum of a superconductor, and has a direct bearing 
on many of the experimentally measurable quantities. 

As an example, consider the standard BCS s-wave su- 
perconductor Q, which has an isotropic gap function 
which is even across the Fermi surface (FS), non-zero 
only in a thin shell of width ~ Tiujd around cf- The en- 
ergy gap is related to the gap function (1) through 
the relation 
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and obeys the celebrated gap equation 
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where 



A^ 2 is the quasiparticle energy, e k is 



the single-particle energy with respect to ep, and Vkk' is 



the effective electron-electron interaction matrix element 
in the reduced BCS hamiltonian. The energy gap can be 
directly observed, for example, through superconductor- 
normal metal or Giaver tunneling. 

To understand the significance of the gap function (Q) , 
let us examine the behaviour of the BCS wavefunction 
close to the FS: 
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which can be rewritten as 
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where the product is now over only half of the k-space 
(e.g. k x > and all k v in 2D). Here we have 



= J_ fc t 

- k ~ V2 [ k 



kt c -ki 



c ki c 



f -k T ) 



(6) 



as the singlet Cooper pair creation operator at points k 
and — k in k-space. It is clear that u£ is the probabil- 
ity amplitude of finding no singlet pair with momentum 
(k, — k), M k w k is that of finding one singlet pair (of charge 
2e), and that of finding two singlet pairs (of total 
charge 4e). The BCS state has identical phase relations 
for various configurations of pair occupancy in k-space. 
That is, when the product in Eq. (||) is expanded out, the 
resulting sum has identical phase for all terms, each term 
corresponding to different configurations of the (k, — k) 
occupancy. Superconductivity can thus be thought of as 
a coherent charge-2e fluctuating state in k-space. Since 
A k is non-zero only in a thin energy shell around the FS, 
the coherent 2e charge fluctuation is concentrated around 
the FS. (Note that this coherence in k-space results in 
phase coherence among the Cooper pairs in real space 
also). Away from the shell, we either have a completely 
filled band (inside the FS) or a completely empty band 
(outside the FS) and hence no charge fluctuations. The 
gap function (]l|), for the BCS wavefunction (||), reduces 
to 



A k = u k w k , 



and is thus a measure of coherent charge fluctuations in 
k-space. 
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In the context of the high-T c cuprates, it is generally 
believed that the gap function is not of simple BCS type 
(s-wave and even across kp). Some of the recent experi- 
ments have suggested that the gap is highly anisotropic, 
possibly with the d x 2_ y 2 symmetry. ]To[ | Yet another pro- 
posal is that of odd-gap pairing, which is a favoured mode 
of pairing in the presence of a strong, purely repulsive 
interaction. || A recent proposal in the context of the 
interlaycr tunneling mechanism is that of anisotropic s- 
wave pairing, which has deep minima but no actual nodes 
on the FS. [^0| However, in spite of intense activity, both 
on theoretical and experimental sides, the issue of the 
symmetry of the gap function still remains unresolved. 
The importance of this issue cannot be overestimated, 
since the gap function is a quantity which, on the one 
hand, has direct experimental consequences, and on the 
other, may offer a unique signature pointing to the right 
kind of mechanism for high-T c superconductivity. 

The aim of the present work is to investigate the be- 
haviour of Ak using exact diagonalization, especially 
with reference to its nodal structure. We would like to 
find out the most likely pairing symmetries, both angu- 
lar and radial, for clusters consisting of single or cou- 
pled chains and planes. In particular, we are interested 
in the d x 2_ y 2 symmetry and odd-gap pairing. In what 
follows, we first present a brief summary of the various 
proposed gap functions, especially d x 2_ y 2 and odd-gap 
pairing (section II). We then present our numerical pro- 
cedure in section III and the results in section IV. Finally, 
conclusions are presented in section V. 

II. D X 2_ Y 2, ODD-GAP AND OTHER PAIRING 
SYMMETRIES 

In order to see, on general grounds, why gap func- 
tions with zeros are preferable as candidates for pairing 
in high-T c superconductors, let us consider the one-band 
large-?/ Hubbard model or the t—J model. These models 
are believed to contain the essential low-energy physics of 
high-T c superconductivity. Let us begin by noting that 
in the large-t/ Hubbard model, the on-site pairing ampli- 
tude An = (cLcJ, ) is diminished considerably on account 
of the restriction on double occupancy. In the t— J model, 
of course, double occupancy is completely projected out. 
Thus Ajj = 0, which reduces, in k-space, to 

A»=^(4 T cL ki )=^A u = 0, (8) 

k k 

where we have used the pairing condition (c^c^,^) = 
Ak<5k.-k'- This is a constraint on the gap function, aris- 
ing from a strong real-space repulsion, and implies the 
presence of zeros in the gap function. 

The constraint, however, is a global constraint, and 
can be satisfied in several ways. One possibility is that 
of odd-gap pairing, in which case the gap function van- 
ishes on the entire Fermi surface, and changes sign across 



the FS. Yet another possibility is that of a simple d x 2_ y 2 
superconductor on a 2D square lattice, for which the gap 
function vanishes along the lines k x = ±k y which in- 
tersect the FS at exactly four points. We give below a 
brief account of these two and some of the other proposed 
pairing states. 

A. Odd-Gap Pairing 

The possibility of odd-gap pairing, in the presence of a 
dominantly repulsive interaction, was anticipated by Co- 
hen H many years ago. The case for odd-gap pairing, 
in the context of high-T c superconductivity, was put for- 
ward by Mila and Abrahams (^] , within the framework of 
the weak coupling BCS theory. They have shown that an 
even solution to the BCS gap equation (^) exists only so 
long as the effective electron-electron interaction is dom- 
inantly attractive, albeit weak. This is the case with the 
standard BCS model of conventional superconductivity, 
where the weak attraction arises due to electron-phonon 
coupling. In the presence of a repulsive interaction, how- 
ever, such a solution ceases to exist when the repulsive 
part of the interaction starts dominating, and an unusual 
solution is shown to exist which is odd across the FS. For 
example, in the toy models considered in as well as 
the energy gap is proportional to (e — ep), which is 
clearly odd across the FS. Moreover, an odd gap is inde- 
pendent of the magnitude of the repulsion, provided that 
it is the repulsion which dominates. This condition is 
met in the large-?/ Hubbard model, and odd-gap pairing 
is thus a candidate for pairing in the high-T c supercon- 
ductors. Mila and Abrahams have further shown that 
certain key features of the tunneling density of states in 
the cuprates, such as the temperature-insensitive peak, 
can be naturally explained on the basis of odd-gap pair- 
ing. This feature is in contrast with s-wave as well as 
d x 2_ y 2 pairing, in which the peak shifts towards ep as 
T^T C . 

Anderson Q has also pointed out in the context of 
the RVB theory that the gap function vanishes on, and 
changes sign across, the FS. For example, in the RVB 
mean field theory at half filling, || the gap function is 
constant and changes sign across the pseudo FS defined 
by cosk x + cosfcy = and satisfies the constraint (||). 
For the doped case, however, the analysis is difficult and 
Anderson did not provide any argument as to why the 
gap should vanish on the FS away from half-filling. || 
In the latter case, it is argued that strong correlations in 
real space lead to a suppression of coherent charge fluc- 
tuations close to the FS, p7| in the sense of Eq. (||), as 
follows. Strongly correlated electrons in ID, described by 
a large-f/ Hubbard model, have certain unique features 
close to the FS. It is well known that in the ID Hub- 
bard model there is singular forward scattering between 
two electrons with opposite spins close to the FS. This 
leads to a finite phase shift [|| at the FS and the conse- 
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quent failure of the Fermi liquid theory, resulting in the 
vanishing of the discontinuity in nk at the FS (Luttinger 
liquid behaviour). It also implies an effective hard-core 
repulsive pseudopotential between electrons with oppo- 
site spins close to the FS. Thus no two electrons close to 
the FS, with opposite spins, can have the same momen- 
tum, thereby making k-points close to the FS essentially 
singly occupied. Coherent pair fluctuations are thus un- 
likely to develop on or very close to the FS, but are not 
forbidden away from the FS. It is therefore likely that 
the line of zeros of Ak implied by the global constraint 
(Eq. (§)) will coincide with the FS. 

Whether these arguments carry over to 2D or not is 
still unclear. It is believed that if the constraint of no 
double occupancy in real space leads to singular forward 
scattering and the consequent failure of Fermi liquid the- 
ory in 2D (in the spirit of Anderson's tomographic Lut- 
tinger liquid picture [^)), a gap function which vanishes 
on the entire FS is a natural candidate for pairing. 

Odd-gap pairing has so far found very feeble experi- 
mental support. However, a very recent work shows 
ARPES data in Bi2212 cuprates in which the gap does 
not vanish along k x = k y , but vanishes on either side of 
this line. This can be interpreted as a possible indication 
that the nodal lines in the gap coincide with the FS. 



B. d x 2_ y 2 Pairing 

The d x 2_ y 2 paired state is characterised by a gap func- 
tion of the form 



Ak = Ad (cos k x — cos k y ). 



(9) 



Support to d x 2_ y 2 pairing comes from numerous sources. 
Theoretical analysis of motion of holes in an antiferro- 
magnetic background shows that exchange of spin fluc- 
tuations can induce singlet pairing of d x 2_ y 2 symmetry, 
for which hole-hole interaction is attractive. || Although 
antifcrromagnetic fluctuations are rather subdued in the 
doped case in comparison with half-filling, it is reason- 
able to expect d x 2_ y 2 pairing as a possibility. 

On the experimental side, numerous experiments have 
been performed recently, which seem to support d x 2_ y 2 
pairing. London penetration depth data from Hardy et 
al. jm shows a low-temperature linear temperature de- 
pendence which is in agreement with a d x 2_ y 2 supercon- 
ductor, and very much in contrast with the BCS s-wave 
case (which shows an exponential behaviour). An exper- 
imental result strongly supporting d x 2_ y 2 pairing is the 
T 3 variation of Cu NMR relaxation rate @. ARPES 
measurements on Bi2212 compunds by Shen et al. fl0(| in- 
dicate lines of zeros compatible with a d x 2_ y 2 state, and 
incompatible with an extended-s state, although these 
measurements are unable to rule out pairings such as 
s + id or an anisotropic s-wave (next section). Certain 
other experiments, which support d x 2_ y 2 pairing, but 
cannot rule out the possibility of other kinds of pairing 



include SQUID experiments on YBCO and electronic 
Raman scattering measurements on Bi2212 compounds. 



On the numerical side, studies in favour d x 2_ y 2 pairing 
include an exact diagonalization study of possible pairing 
symmetries by Riera and Young []l5f (on which we com- 
ment in the next section). Dagotto and Riera jl6| find 
bound hole pairs with d x 2_ y 2 symmetry near half-filling 
and on the verge of phase separation, although the value 
of J/t ~ 3 is rather too large in comparison with the 
experimentally relevant parameter range J/t ~ 0.3. An- 
other notable work calculates an anomalous Green's 
function for one and two holes using Lanczos technique, 
and finds signals of d x 2 __ y 2 pairing over a wide range of 
parameter J/t (not specified). 



C. Other Proposed Pairing Symmetries 

A variational Monte Carlo study by Li et al. (Hf l 
supports s + id symmetry rather than a pure d x 2_ y 2. 
This calculation compares RVB variational states with 
extended-s, d x 2_ y 2, and mixed states s + d and s + id 
symmetries. The authors show that the mixed states are 
energetically preferred over pure s or d x 2_ y 2 states, and 
that the s + id state overcomes finite-size effects at a lat- 
tice size smaller than that required by the d x 2_ y 2 state. 
The authors also argue that the Knight shift data in the 
YBCO cuprates is better explained by s + id rather than 
pure d x 2_ y 2 pairing. Wheatley |l9| has also analyzed the 
stability of various mixed states involving d x 2_ y 2. 

Chakravarty et al. |^0| have recently proposed an 
anisotropic s-wave pairing state based on the interlayer 
tunneling mechanism [|T) . This state has no actual nodes 
on the FS although deep minima occur, and the gap is 
highly anisotropic. Since ARPES results, such as |[(|, are 
unable to detect the sign of the gap function (as well as 
the presence of a node unambiguously owing to the cur- 
rent energy resolution), they are consistent with such an 
anisotropic s-wave pairing as well. Interlayer tunneling is 
also responsible for a temperature-dependent anisotropy 
in the gap, as shown by Muthukumar and Sardar p3| , 
which is incompatible with the d x 2_ y 2 scenario. Such a 
temperature-dependent anisotropy has been observed in 
the ARPES data by Ma et al. pi 



III. NUMERICAL PROCEDURE 

In our exact diagonalization computation, since we 
have chosen to work in a number-conserving basis, it is 
not possible to evaluate Ak directly. 24 We thus define 
a correlation function Ak.k' as 



(10) 



where the average (. . .) is the expectation value in the 
ground-state of a finite cluster, obtained via exact diago- 
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nalization. This quantity is clearly a measure of coherent 
pair fluctuations between states (k, — k) and (k', — k') in 
the ground state. 

We find the gap function A k by diagonalizing this two 
particle reduced density matrix ( |f0| ) , which has the eigen- 
function decomposition 

A k)k > = ^A a A Q (k)A*(k'). (11) 

a 

Here, X a and Ao,(k) are the a th eigenvalue and eigen- 
function of the N x N matrix . Here N is the num- 
ber of points in the Brillouin zone and the index a orders 
A a 's as Ai > A2 > A 3 > . . .. Superconducting ODLRO is 
signalled J25| by a macroscopic separation of the largest 
eigen value Ai from the next one A2, i.e. Ai — A2 ~ N. 
The required gap function Ak is the eigenvector Ai(k) 
corresponding the largest eigen value Ai, which is the 
Cooper pair state with least energy for which condensa- 
tion will take place, whereas the other eigenvectors rep- 
resent the various excited states of the Cooper pair. For 
example, for the standard BCS ground state, Ai = N , 
A x (k) = Mk^k and X a = 0, for a = 2, 3 . . . N. 

A similar procedure was adopted by Riera and Young 
jl5| in the context of the t — J model with the three- 
site term. However, the authors argue that from the 
expression (|l0|) translated into real space, 

A ** = ^2 E e !k - (Rfc - R!) e - 4k '- (R '"- R " ) <&Lw, 

klmn 

(12) 

(where tf mn = ^(c^c^ -c^4 T ) creates a spin-singlet 
pair of electrons at sites m,n), one should retain only 
those terms which have all the four indices k, I, m, n dis- 
tinct, since it is only these terms which correspond to true 
pairing terms, rather than mere charge and spin corre- 
lations. The authors post facto justify this procedure on 
the basis that the d x 2_ y 2 symmetry starts dominating 
very close to where the two-hole binding energy becomes 
negative, as a function of J. We do not agree with this 
argument, since there is no clear cut way to justify such 
a removal a priori, in our calculation, since we are calcu- 
lating a quantity in k-space; Eq. (|l2|) then dictates that 
all terms be retained in the Fourier transform. Indeed, 
such removal in the U < case leads to an incorrect 
symmetry, whereas the gap function is known to be of 
a simple BCS type, i. e. an even function across the 
FS (with an s-wave symmetry in 2D). This is illustrated 
in Fig. 1 in which we plot the gap function correspond- 
ing to the highest eigenvalue for an 8-site Hubbard chain 
with U = — 5, with and without following this removal 
of terms. Without the removal of terms, we see the ex- 
pected behaviour, whereas with the removal, we observe 
a sign change across ki?. On a y/8 x y/8 Hubbard plane 
at U = —5, with the removal of terms, we have observed 
a double sign change across the FS, indicating two nodal 



lines encircling the FS on either side, whereas without 
such removal, we see the expected BCS-like symmetry. 

We shall thus present results without the removal of 
these "on-site" terms. Note that with removal, we indeed 
obtain similar results as Riera and Young. |l9] We would 
also like to point out that unlike their work, the focus 
of our work is the nodes in the gap function. We shall 
refer to the gap function corresponding to the highest 
eigenvalue of the density matrix ([l0]) as the "topmost" 
gap function, and the term eigenspectrum will be used to 
refer to the eigenvalue spectrum of the density matrix. 

IV. RESULTS AND DISCUSSION 

The clusters used for these calculations are single-layer 
clusters of 8 (Hubbard) and 16 (t — J) sites, and clusters 
with two 8-site t — J layers coupled by an interlayer term 
of the form 

H ± = -t± E c 1 - "U)4c(i - <_ a ) (13) 

ilraa 

where l,m are layer indices (= 1,2), i is the site index 
within a layer and a is the spin label. 

A. Results in One Dimension 

We first present the gap function for the ID clusters. 
Fig. 2 displays the topmost gap function for U — 10 at 
various fillings. A clear change of sign is visible across 
the corresponding Fermi point (§ at half filling, j for 2 
and 4 holes), indicating the presence of a node near kp . 
Similar behaviour is observed for all U > 0, and over a 
range of fillings (0, 2, 4 and 6 holes). 

The results for a 16-site t—J chain are presented in Fig. 
3 and 4, and are representative of the parameter range 
J = 0.08 — 0.32. Fig. 3 presents eigenspectrum of the 
density matrix for various cases: half-filling at J = 0.24 
(1), 2 holes at J = 0.32, 0.24, 0.16, 0.08 (2)-(5). All these 
spectra display a prominant, well-separated largest eigen- 
value, indicating the possibility of ODLRO (as far as a fi- 
nite size calculation can reveal). For comparison, we have 
also plotted the eigenspectrum for a antiparallcl triplet 
Cooper pair density matrix in Fig. 3 (6) for J — 0.24, 2 
holes, which does not show a well-separated eigenvalue. 
Fig. 4 displays the topmost gap function at half-filling 
and 2 holes, where a clear sign change across k^ indicates 
a node near k^?. 

Finally, Fig. 5 plots the topmost gap function for a 
single 8-site chain in a system of two 8-site t — J chains 
coupled by an interlayer coupling term (|l3|), for 2 holes 
with J = 0.24 and t± = 0.05 — 0.7. We see again a clear 
node near ki? = J very similar to the case of a single 
16-site t — J chain. 
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B. Results in Two Dimensions 

We now turn to the more interesting case of 2D clus- 
ters. These results are obtained for V° x v° Hubbard, 
4 x 4 t - J and coupled \/8 x \/8 t — J planes. The 
results are representative of a range of parameter val- 
ues and fillings as specified below. V8 x y/8 Hubbard 
cluster: all U > (0, 2 and 4 holes); 4 x 4 f - J plane: 
J = 0.08, 0.16, 0.24, 0.32 (0 and 2 holes); coupled V8 x v 7 ! 
f - J planes: J = 0.24,0.32, i ± = 0.05 - 0.7 (0 and 2 
holes). To bring about the symmetry features of a gap 
function, we shall display the numerical values of the co- 
efficients of the corresponding eigenvector (of the density 
matrix (^0[)) at all k-points in the Brillouin zone. The FS 
(and the nodal lines of d x 2_ y 2, wherever required) will be 
indicated by dashed lines. 

We first discuss the characteristics of the eigenspec- 
tra. In all the cases presented, the eigenspectrum dis- 
plays a separated largest eigenvalue, although not as well- 
separated as in the ID case. Fig. 6 (l)-(3) are eigenspec- 
tra for 2 holes, respectively for a y/8 x y/8 Hubbard cluster 
with U = 10, 4 x 4 i — J plane with J = 0.24 and coupled 
\/8 x y/8 t — J planes with J = 0.24, t± = 0.2. The most 
striking feature of these results is the close interplay be- 
tween the odd-pairing and the d x 2_ y 2 gap functions. The 
odd-paired state is the topmost for all the three clusters, 
with the next dominant state as the d x 2_ y 2 state on 4 x 4 
and coupled \/8 x \/E clusters. For the y/8 x y/8 Hub- 
bard cluster with 2 holes, however, the d x 2_ y 2 state is 
quite low down in the eigenspectrum. These two states, 
for a 4 x 4 plane, are displayed in Fig. 7. 

C. Discussion 

The overall picture that emerges out of this study is 
as follows. In ID and 2D, the gap function shows the 
presence of a node near or at kp , indicating an odd- 
paired behaviour. The d x 2_ y 2 state is often seen as the 
next dominant state. Quite clearly, the numerical results 
seem to indicate the presence of a nodal "surface" (GNS) 
in the gap function closely following the FS. Recall that 
in the RVB mean-field theory at half-filling, the GNS 
coincides with the FS in 2D (also seen in our numerical 
results (not shown)). It may be argued that as doping 
increases from zero, the FS shrinks, dragging part of the 
GNS with it. |27| The non-overlap of the GNS with the 
FS in some parts in the BZ is possibly caused due to 
the enhancement of the interlayer pair tunneling matrix 

element t± ^ in those directions in k-space. For example, 

as emphasized by Chakravarty et al. |Q , tkitl [ s largest 
in the (0, tt) and (tt, 0) directions which enhances pairing 
in those regions of the FS by keeping the GNS away. Fig. 
8 incorporates this idea as well as the eight nodal points 
of the recent ARPES data 0. Based on the present 
numerical work, we believe that it may be possible to 



construct novel scenarios where such a nodal surface is 
also consistent with d x 2_ y 2 pairing. This and such other 
aspects are discussed elsewhere. [E7| 

The major handicap of any exact diagonalization cal- 
culation is the cluster size restriction. It is not possi- 
ble for us to rule out finite size effects, but the features 
we have observed are consistent within the sizes and ge- 
ometries considered. It is indeed possible that finite size 
effects are more serious when symmetry issues are con- 
cerned, in comparison with, for example, the finite size 
effects in simple spin, charge or pairing correlation func- 
tions, u 

Note that we do not interpret the present results as 
indicative of the presence of superconducting ODLRO 
in purely 2D systems. However, it does indicate that 
these pairing correlations show an odd-paired behaviour 
across the FS, a feature which seems to persist even in 
the presence of an interlayer coupling term. 

V. SUMMARY 

To summarise, we have investigated the behaviour of 
the gap function At across the Fermi surface in the doped 
t — J and large- U Hubbard clusters in one and two di- 
mensions, with special reference to odd-gap and d x 2_ y 2 
pairing. For this purpose, we diagonalize a reduced two- 
particle (Cooper pair) density matrix computed through 
exact diagonalization. The largest eigenvalue of this den- 
sity matrix is well-separated from the rest, indicating the 
possibility of ODLRO. The corresponding eigenvector, 
the Cooper pair wavefunction with least energy, reflects 
the symmetry of the gap function At- The results on 
8- and 16-site single or coupled chains and planes clearly 
indicate a change of sign of Ak across kp, indicating the 
presence of a node close to the entire FS, a signature 
of odd-pairing behaviour. This state, on single and cou- 
pled planes, in addtion, is isotropic, indicating an angular 
symmetry which is of the s-wave type. We also observe a 
close interplay between the odd-paired and d x 2_ y 2 states 
in the 2D systems. 
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Figure Captions 

1. Topmost gap function for an 8-site Hubbard chain 
with 2 holes at U — — 5, with and without the 
removal of terms from the gap expression (p^); 
k F = f . Notice the sign change across the FS 
with the removal of terms. 

2. Topmost gap functions for an 8-site Hubbard chain 
with U = 10, for 0, 2 and 4 holes; k F = § , f , f 
respectively. 

3. Eigenvalue spectra of the density matrix for the 16- 
sitc t - J chain. (1) Half-filled, J = 0.24, (2)-(5) 2 
holes, J = 0.32,0.24,0.16,0.08, (6) triplet antipar- 
ailel eigenspectrum for 2 holes, J = 0.24. 

4. Topmost gap functions for a 16-site t — J chain at 
half filling (J = 0.24) and 2 holes (J = 0.08-0.32). 
k F = § , for and 2 holes respectively. 

5. Topmost gap function for a single 8-site t — J chain 
in a coupled system of two chains with 2 holes, for 
values of t± = 0.05 - 0.7, J = 0.24. k F = § for 

ky = 0. 

6. Eigenvalue spectra of the density matrix for the 
three 2D clusters at half filling and for 2 holes. (1) 
s/& x V8 Hubbard plane, 2 holes, U = 10; (2) 4 x 4 
t - J plane, 2 holes, J = 0.24; (3) coupled V& x y/E 
t - J planes, 2 holes, J = 0.24, t ± = 0.05 - 0.7. The 
labels "o" and "d" near a level respectively stand 
for odd-paired and d x 2_ y 2 states. A numeral near 
a level indicates the degeneracy. 

7. The topmost two degenerate gap functions for a 
4x4 plane with 2 holes, J = 0.24. Points marked 
by squares lie within one Brilliouin zone, diamonds 
belong to the neighbouring ones. The FS and the 
nodes of the d x 2_ y 2 state are indicated by dashed 
lines. The numbers above and below a square are 
the coefficients of the gap function respectively for 
the topmost (odd-paired) and next (d x 2_ y 2) states. 

8. A suggested form of the GNS with s-symmetry. 
The solid line is the FS and the dotted line is the 
GNS. The + and - symbols indicate relative signs 
of the gap function across the GNS. 
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